Single scatterings in single artificial atoms: Quantum coherence and entanglement 
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We employ the quantum-jump approach to study single scatterings in single semiconductor quan- 
tum dots. Two prototypical situations are investigated. First, we analyze two-photon emissions 
from the cascade biexciton decay of a dot where the single-exciton states exhibit a fine-structure 
splitting. We show that this splitting results for appropriately chosen polarization niters in an oscil- 
latory behavior of two-photon correlations, and carefully examine the proper theoretical description 
of the underlying scattering processes. Secondly, we analyze the decay of a single-electron charged 
exciton in a quantum dot embedded in a field effect structure. We show how the quantum properties 
of the charged exciton are transferred through tunneling and relaxation to the spin entanglement 
between electrons in the dot and contact, and identify the pertinent disentanglement mechanisms. 
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I. INTRODUCTION 

Ultrafast semiconductor spectroscopy provides a 
unique laboratory for the investigation of coherence and 
scattering effects in solids: 1,2 a short laser pulse creates 
electron-hole pairs propagating in phase with the excit- 
ing laser; this initial phase coherence is consecutively de- 
stroyed through mutual interactions among the photoex- 
cited carriers and couplings to environmental degrees of 
freedom, e.g. phonons, thus leading to decoherence 3,4 
which can be monitored again by means of ultrafast spec- 
troscopy. Recent years have seen spectacular examples of 
such ultrafast semiconductor spectroscopy, e.g. the ob- 
servation of phonon quantum beats 2 ' 5 ' 6 or the buildup of 
screening, 7 and have revealed dephasing and relaxation 
times of the order of femto- and picoseconds for carrier- 
carrier and carrier-phonon interactions, respectively. 

In semiconductor quantum dots, 8,9 sometimes referred 
to as artificial atoms, the strong confinement in all three 
spatial directions results in a substantial suppression of 
decoherence as compared to higher-dimensional semicon- 
ductors, i.e., bulk, quantum- wells, or quantum wires. On 
the one hand, Coulomb interactions among carriers cap- 
tured in the dot do not result in scattering but only give 
rise to energy renormalizations of the electron-hole few- 
particle states. 10-15 On the other hand, at low temper- 
atures phonon-mediated dephasing of the electron-hole 
states of lowest energy is of only minor importance, 16 ' 17 
and can become even negligible in comparison to sponta- 
neous emission of photons which occurs on a nanosecond 
time scale. 18 ' 19 These remarkable features render quan- 
tum dots ideal candidates for the solid-state implemen- 
tation of quantum engineering, 20,21 e.g., for the purpose 
of quantum computation, which has led to considerable 
research activities in recent years. 22-30 

Observation of optical coherence effects in ensembles of 
dots is usually spoiled by inhomogeneous line broaden- 
ing due to dot size fluctuations, with typical broadenings 
comparable to the level splittings themselves. To over- 
come this problem, within the last couple of years a num- 



ber of experimental techniques were developed to allow 
the observation of single quantum dots, hereby establish- 
ing the rapidly growing field of single-dot spectroscopy. 31 
Such single-dot measurements are not only indispensable 
for the observation of quantum coherence and the im- 
plementation of quantum engineering, but also open the 
possibility to study decoherence and relaxation in sin- 
gle quantum systems. The theoretical analysis of sin- 
gle quantum system, however, is more cumbersome as 
compared to ensembles of corresponding systems. In the 
latter case, one usually employs the framework of quan- 
tum transport theory or quantum kinetics 1,2,32 which de- 
scribes the time evolution of macroscopic quantities, e.g., 
the interband polarization, directly accessible to experi- 
ment; here, the detailed evolution of one specific subsys- 
tem is completely irrelevant, which allows to submit the 
problem to the laws of statistical physics. On the other 
hand, for a single system one has to be more specific 
about the occurrence of single scatterings and their mon- 
itoring through a measurement apparatus. The ques- 
tion of how to theoretically describe such problems first 
arose almost two decades ago when it became possible to 
store single ions in a Paul trap and to continuously mon- 
itor their resonance fluorescence, and led to the devel- 
opment of the celebrated quantum-jump approach. 33 ^ 35 
This approach combines the usual master-equation ap- 
proach with the rules of demolition quantum measure- 
ments, 35,36 and provides a flexible tool for the description 
of single-system dynamics subject to continuous monitor- 
ing. 

In this paper we employ the quantum-jump approach 
to two representative examples for the observation of sin- 
gle scatterings in single semiconductor quantum dots. 
First, we study the cascade decay of a biexciton, which 
has recently attracted considerable interest: 37-40 in a first 
step, the biexciton decays through emission of a pho- 
ton to the spin-degenerate exciton states; this final-state 
ambiguity leads to an entanglement of photon and exci- 
ton; when subsequently the exciton decays, the photon- 
exciton entanglement is transferred to a two-photon en- 
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tanglement. 37 ' 39 As shown below, for dots with an exci- 
ton fine-structure splitting 41 the two-photon correlation 
function of this cascade decay exhibits an oscillatory be- 
havior, provided that the photons are detected in a prop- 
erly chosen polarization basis (see Ref. 42 for a similar 
experiment). Our analysis carefully examines the pho- 
ton emission processes, and improves upon the simple- 
minded framework of Fermi's golden rule to provide a 
proper description of the problem under concern. Our 
second example is based on an experiment recently per- 
formed by Zrenner et al. 43 in which the authors used 
a quantum dot embedded in a field-effect structure to 
transform an optically excited exciton through tunnel- 
ing into a photocurrent. Here, we study the decay of 
a charged exciton, i.e., of a Coulomb- renormalized com- 
plex consisting of two electrons and a hole. 13-15 Through 
the tunneling process one electron with either spin-up 
or spin-down orientation is transferred from the dot to 
the reservoir; since the electron in the dot has an op- 
posite spin orientation, the spins of the reservoir- and 
dot-electron become entangled. 44 In contrast to the op- 
tical biexciton decay additional difficulties arise because 
the tunnel-generated electron and hole do not propagate 
freely, as photons would in the corresponding scheme, but 
are subject to interactions in the contact. We develop 
a prototypical, though simplified, description scheme of 
the scattering channels in the reservoir, and demon- 
strate that this spin entanglement is robust against spin- 
unselective scatterings and is only deteriorated by spin- 
selective scatterings. 

We have organized our paper as follows. In Sec. II we 
discuss the monitored biexciton cascade decay in a sin- 
gle quantum dot. We specify our model system, provide 
a short discussion of the master-equation and quantum- 
jump approach, and finally show how to obtain the sys- 
tem's time evolution through unraveling of the master 
equation. The charged-exciton decay of a dot embedded 
in a field-effect structure is discussed in Sec. III. We de- 
velop a theory accounting for the complete cascade pro- 
cess of the buildup of coherence through tunneling, the 
swapping of quantum coherence to the spin entanglement 
through dephasing and relaxation, and the process of dis- 
entanglement through spin-selective scatterings. Finally, 
in Sec. IV we draw some conclusions. 



II. TWO-PHOTON CORRELATIONS 
A. Quantum dot level scheme 

1. Linear basis 

We examine a prototypical level scheme depicted in 
Fig. lb, which consists of the groundstate (no electron- 
hole pairs present), the spin-degenerate single-exciton 
states s = 1,2 of lowest energy, and the biexciton state 
3; other dot states, such as excited excitons or biexci- 
tons, are energetically well separated and can be safely 
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FIG. 1: (a) Schematic picture of the proposed setup. Photons 
from the biexciton cascade decay are monitored by the pho- 
todetectors Q s and Q , where Q s (fio) is sensitive to photons 
from the 3-s (s-0) decay, (b) Level scheme as described in the 
text: is the groundstate, s = 1, 2 the spin-degenerate single- 
exciton states, and 3 the biexciton one. The arrows indicate 
the optically allowed transitions and the symbols indicate x- 
or y-polarization for photons propagating along z. 



neglected. The single-exciton energies e s = E a ± 5/2 are 
assumed to be split by a small amount 6 (typically a few 
tens of /xeV) because of the electron-hole exchange in- 
teraction 41 and the corresponding states to be polarized 
along x and y. Finally, the biexciton energy €3 = 2E a — A 
is reduced by the biexciton binding energy (typically a 
few meV) because of Coulomb correlations. 12 Within 
this scheme, the Hamiltonian of the quantum-dot states 
reads: 



H d Q = e 5 | S )( S | + e 3 |3)<3|, (1) 

s=l,2 

where energy zero is given by the groundstate. In what 
follows, we shall make use of the fact that the groundstate 
is optically connected to the biexciton state 3 through 
two independent pathways, namely 0-1-3 and 0-2-3, 
which allows to establish a quantum coherence between 
the single-exciton states directly observable in optical ex- 
periments. 

Because of the important role of states s = 1, 2 we shall 
find it convenient to introduce the Pauli matrices a± = 
|1)(2| + |2)(1|, a 2 = -i (|1)(2|-|2)(1|), a 3 = |1)<1|-|2)<2|, 
which together with the unit matrix 1 = |1)(1| + |2)(2| 
form a complete basis within the single-exciton subspace. 
The Hamiltonian of Eq. (1) can then be rewritten as 



H« = E l + ~a 3 +e 3 \3)(3\. (2) 
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2. Circular basis 

In the linear basis the Hamiltonian (2) is diagonal. 
Hence, a laser pulse propagating in z-direction with po- 
larization along x or y allows to selectively excite one of 
the two single-exciton eigenstates. On the other hand, 
any other light polarization will induce a superposition 
state, e.g., 

|+) = (|l)+z|2»A/2 

I-) = i(|l)-i|2»/>/2 (3) 

for circular polarization a± . Further below we shall find 
it convenient to switch between the linear basis s = 1, 2 
and the circular one, Eq. (3), through a unitary transfor- 
mation {/=(! + ioi)/\/2, where 

U(E l + l(T 3 )Ui =E l+ 5 -a 2 . (4) 

In other words, the energy splitting 5 translates in the 
circular basis to a coupling between states a = ±. Thus, 
if the system is initially prepared in, e.g., state |+) (which 
is not an eigenstate of the Hamiltonian), it will start to 
oscillate between \a — ±) where the oscillation period is 
given by 1/8. 

B. Light-matter coupling 

We describe the light-matter coupling within the usual 
dipolc and rotating-wave approximations 45,46 

#op = T + ■ (5) 

Here, 6^ is the electric field of the light evolving with 
positive or negative frequency components, 45 ' 47 

V = ^e s (\s)(0\ + \3)(s\) (6) 

s 

the interband polarization of the quantum dot, the 
dipole moment, e s a unit vector along x or y, and s 
denotes the state with polarization perpendicular to s. 

1. Coherent optics 

Suppose that the quantum dot initially in its ground- 
state is excited by a short laser pulse. We assume that 
6^ is sufficiently intense to be treated classically, and 
that the spectral width of the pulse is sufficiently broad 
to excite both 0-s transitions with comparable strength 
but narrow enough to inhibit transitions between s and 
3. Then, such a pulse with circular polarization a + will 



induce a superposition — a\0) + /3|+), where a and 
f3 are determined by the laser parameters. Since |+) is 
not an eigenstate but a coherent superposition of states 
s = 1,2, as time goes on the system will oscillate between 

l±>- . , 

This optically induced quantum coherence corresponds 

to an oscillating interband polarization which can be di- 
rectly monitored by a second, time-delayed laser pulse 
(see, e.g., Refs. 19,48 for the observation of such coher- 
ence beats). Let us assume that the frequency of the 
second pulse is tuned to the s-3 transitions and that the 
delay time is short as compared to the decoherence times. 
Then, the absorption of the second pulse will sensitively 
depend on the time delay between the two pulses, and 
will exhibit oscillations on a time scale of 1/8. 



2. Incoherent optics 

In this paper we shall address the question whether 
it is possible to observe a similar effect by means of in- 
coherent optics. We assume that initially the dot is in 
the biexciton state 51 and consider the setup depicted in 
Fig. la, where the optical decay of photons along ±z is 
monitored by two photo detectors. For this geometry the 
usual optical selection rules (to be discussed below) ap- 
ply. In addition, we assume that photo detector f2 s (f2 G ) 
is sensitive to photons from the 3-s (s-0) decay, 38 which 
can be achieved through spectral filtering (and restricts 
the time resolution to ~ 1/A <C 1/8) . Suppose that at 
time t a photon with polarization cr_ is detected from the 
biexciton decay. Analogously to the coherent excitation- 
scenario described above, we might assume that after 
photon detection the system is prepared in state |+). 
Since this is not an eigenstate of the Hamiltonian, the 
system will oscillate between |±) which could be detected 
by monitoring the probability distribution for the second 
photon emission. Thus, time-resolved photon correlation 
measurements together with appropriate frequency and 
polarization filtering should provide a means for the mea- 
surement of quantum coherence in an incoherent optical 
experiment. 

At this point we encounter a conceptual problem. If 
we would describe the photon emissions through Fermi's 
golden rule we would obtain scattering rates T 3 ^\ be- 
tween the asymptotic states before, 3, and after, A, the 
scattering. Here, the photon correlations would strongly 
depend on the choice of the basis A: in the circular basis, 
the detection of a ct_ photon would correspond to an opti- 
cal decay from 3 to |+) and would result in the aforemen- 
tioned coherence beats. On the other hand, in the linear 
basis the detection of a er_ photon had to be associated 
with equal probability to either the 3-1 or 3-2 decay, and 
the state after photoemission would be a mixture with- 
out any quantum coherence between s = 1,2 (since the 
nature of a scattering is intrinsically incoherent). Here is 
the problem: within the framework of Fermi's golden rule 
different predictions arise for different basis sets. This is 



4 



a very unfortunate situation which calls for a more care- 
ful theoretical analysis. It is the purpose of the remainder 
of this section to reexamine the problem of photon cor- 
relation measurements from the biexciton decay and to 
provide a comprehensible answer. 

3. Photon-matter coupling 

We describe the coupling between the quantum-dot 
states and photons through Eq. (5), where the electric 
field of a photon with wavevector k and polarization a 

(+) . / 27TW fc \ 5 „ 

£ ka ■ efeffflfeCT (7) 

is expressed in terms of the usual bosonic field operators 
afcer- 45 ' 47 ' 52 Here, Lo k = ck/n is the light frequency, c the 
speed of light, k the semiconductor dielectric constant, 
n = y/n its refractive index, and e ka the polarization 
vector (see also appendix A). With Eq. (7) the light- 
matter coupling (5) can be cast to 

H o P = * E (.^-«L (|0><*| + |5><3|) - h.c.) , (8) 

fc,(TS 

with the coupling constant gk,as = (2nuJk / n) 5 /j, {e* ka e s ). 
Additionally, in the following we will need the Hamilto- 
nian for free photons H2 = J2ka a\. a a k(T . 

C. Environment coupling 

Next, we address the problem of how to describe the 
time evolution of a single quantum system in contact with 
its environment, i.e., an open system.' 15,53 ' Suppose that 
the Hamiltonian describing the system coupled to its en- 
vironment can be split into the parts H a + V, where H Q 
accounts for the uncoupled system and reservoir (here 
H£ + H2) and V to the coupling (here -ffj p ). Because of 
the incoherent nature of such system-environment cou- 
plings we have to adopt a density-matrix description, 
where the density operator p fully characterizes the sys- 
tem: its diagonal elements provide information about 
state occupancies whereas the off-diagonal elements ac- 
count for quantum coherence. In lowest-order time- 
dependent perturbation theory and assuming that at 
time zero system and environment are effectively decou- 
pled, the time evolution of p is governed by: 45 ' 53 ' 54 

ft = fdt'lv'mv'inpt]], (9) 

Jo 

where the superscript I denotes the interaction repre- 
sentation according to H Q . Two points are worth em- 
phasizing: first, Eq. (9) is of Markovian form, i.e., the 



time evolution of p at time t is completely determined 
by the density operator itself; as discussed in detail in 
Ref. 54, this approximation is appropriate for systems 
where higher-order contributions to Eq. (9) are negligi- 
ble, or in other words, where scatterings occur seldomly 
on a time scale given by H Q (here ~ 1/5); secondly, the 
time integral in Eq. (9) extends to the past of the system 
and thus explicitly accounts for the buildup of scatter- 
ings. Note that the limit t — > oo in Eq. (9), which we 
will not perform below, would recover the framework of 
Fermi's golden rule and would lead to the aforementioned 
final-state ambiguities. 

Below we will need the expression 




X3fc,«4a I S u »-*rf + \s){3\ e *(«*+'»-«3)t'] 

-9t,, s a ka [ \s)(0\ e-^"-^ + \3)(s\ e -*(«*+«.-«.)*'] ) , 

(10) 

which accounts for the buildup of scatterings. The time 
integrals of the exponentials can be performed analyti- 
cally, 

ft Hit _ I 

/ dt'e int = — = e j0t 7 (r!,t), (11) 

Jo « 2 

with 7(fi,t) = sin Sit /CI, and we have neglected in the 
last expression on the right-hand side terms which only 
contribute to an energy rcnormalization but not to de- 
phasing and relaxation. 53 7(0, t) is a symmetric function 
of Cl which in the limit t — > oo gives Dirac's delta func- 
tion. Finally, from Eq. (10) and the definition for 7(0, t) 
we obtain 




xgk,*sa k o [I°)( s l7(^fe -e s ,t)+ \s)(3\j(uj k + e s - e 3 ,t)] 
-9*k,as a ka [\s)(0\j(uj k - e s ,t) + \3)(s\j(uj k + e s - - e 3 ,t)] ) . 

(12) 

This expression provides the starting point for our fol- 
lowing discussion. 

1. Unmonitored decay 

We shall now be more specific about the question what 
happens to the emitted photon. Quite generally, the pho- 
ton can propagate into a direction not covered by the 
photo detectors, in which case it will not be observed, or 
it can be measured by one of the two photo detectors, in 
which case we will acquire additional information about 
the state of the system. 
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Let us first discuss the case of the unmonitored de- 
cay, i.e., photon emission into a segment not covered by 
the photo detector. For simplicity, we only consider the 
biexciton decay monitored by photo detector Q s , where 
analogous conclusions hold for the exciton decay mon- 
itored by £1 . In the spirit of the usual framework for 
the description of open-system dynamics, 45,53 we shall 
trace out the reservoir degrees of freedom. This proce- 
dure implies that we keep ignorant about the state of the 
reservoir. 4 Denoting with 7 g il s the photon modes not 
covered by the photo detector, relaxation and dephasing 
in the dot due to unmonitored decay are described by the 
expression 

pf 1 - -tr 75ZOs /* dt' [V T (t), [V'it'^pt] ] , (13) 
Jo 

where p d denotes the reduced density operator of the 
quantum dot. Eqs. (12,13) give us the prescription of 
how to compute the generalized scattering rates. Since 
each system-reservoir coupling V comes with the bosonic 
field operators a and a) , the resulting expressions contain 
expectation values of aa, aW, aa^ , and a^a. We shall 
treat the reservoir (photons) as a thermal bath at zero 
temperature, i.e., we assume the usual correlation func- 
tions 45 - 53 

tr iP a k<r a k'o> — tr 7/ 9a fc(T a k ' a ' = (14a) 
tr 7 pa fc(T a\, a , ^ 5 kk >5 aa i . (14b) 

With these correlations it becomes possible to decompose 
Eq. (13) into two parts, which, in the spirit of the Boltz- 
mann equation, can be considered as generalized out- 
and in-scatterings. We start by evaluating the double 
commutator [V, [V, p\] = VV'p + pV'V - VpV - V'pV. 
When tracing over the photon degrees of freedom, we re- 
alize that the only non- vanishing terms are proportional 
to aa^p, paa\ and a^ pa. Here, the first two terms do 
not change the number of photons. Correspondingly, we 
shall ascribe the contributions VV' p + pV'V to gener- 
alized out-scatterings which lead to decoherence and re- 
laxation in the dot. If only these terms would be con- 
sidered in the time evolution of p d , the trace of the den- 
sity operator, which describes the probability of finding 
the system in any of the quantum-dot states, would de- 
crease. This lack of norm conservation is restored by 
the remaining terms VpV' + V' pV which describe the ef- 
fects of photon emissions, and may thus be interpreted as 
generalized in-scatterings. As will be shown in Sec. HE, 
this decomposition into out- and in-scatterings allows for 
an efficient solution scheme through an unraveling of the 
master equation. 53 

2. Monitored decay 

The case of a monitored decay has to be treated with 
slightly more care. Quite generally, such monitoring re- 



quires to specify in more detail the quantum mechan- 
ical measurement process. To this end, in the follow- 
ing we shall employ the framework of the quantum-jump 
approach. 35 ' 36 ' 55 We assume that the photon counting 
can be approximately described by a series of repeated 
gedanken measurements at times t n — nA t . Here, A t 
is assumed to be sufficiently short to allow for a time- 
resolved measurement of the system's dynamics, i.e., 
At <C l/<5. On the other hand, A t should be sufficiently 
large to enable the buildup of complete scatterings, such 
that the measurement process itself does not inhibit or 
enhance scatterings (quantum Zeno effect). 35 ' 36 

Suppose that at time t the density operator p t cor- 
responds to an uncoupled dot and reservoir. With the 
projector P^ s associated to the result of the measure- 
ment 3 ' 36 (no photon count, n = 0, or photon count, 
n = 1), the density operator at a later time can be eval- 
uated to 35 - 36 - 55 

Pt — -P" 8 fdt'W'iWVlfH]]^-. (15) 
Jo 

If within the time interval A t no photon is detected from 
the biexciton decay, either no photon has been emitted 
(which corresponds to a projector-like measurement Pq s 
on the photon vacuum) or a photon has been emitted into 
a segment f2 £ fl s . On the other hand, photon detection 
in A t corresponds to a projector-like measurement V^ s on 
the single-photon subspace. Through the measurement 
we know that the biexciton has decayed, and accordingly 
we have to modify the density operator to account for 
this increase of information. 

We can again evaluate the double commutator 
[V, [V',p]] in Eq. (15) and split the resulting expressions 
into generalized out- and in-scatterings; however, within 
the context of monitoring these terms acquire a different 
interpretation: as regarding out-scatterings, anonvanish- 
ing contribution only arises for projection on the photon 
vacuum, which becomes tr 76 n s (VV p+pV'V). This term 
alone is responsible for a decrease of Xvp d associated to 
the diminishing probability of finding the system in the 
initial biexciton state. 35 To restore this lack of norm con- 
servation, we have to additionally consider in-scatterings. 
The corresponding contribution, which is associated to 
photo detection, is non-zero only if projected on ¥ 1 3 , 
and becomes ti J& a s (VpV + VpV). This term gives us 
the prescription of how the density operator must change 
after photon detection. 35 ' 36 



D. Master equation 

With the results of the previous subsections we are 
now in the position to improve upon the simple-minded 
Fermi's golden rule framework. First, we note that for 
both unmonitored and monitored decay the resulting ex- 
pressions are, e.g., of the form ti^^Q'VV' p, where the 
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photon modes extend over a given segment O'. Let us consider one specific contribution 



9*k,*s9k,<Ts>l{uk ~ z s > ,t) = V -n 2 (e ka e s )(e* k<J e s ,)~/(L; k - e s ,,t) 

hi 

ka ka 

2 7J -„2 r roc 

= (2tt) _3 / dn(e ka e s )(e* krT e s <) / k 2 dkuj k j(uj k - e s > ,t) , 

K J new Jo 



(16) 



r 



where we have explicitly evaluated the trace over the pho- 
ton degrees of freedom and used the definition of Eq. (8) 
for the coupling constants g k .as- The quantity of primary 
interest is the integral over k, which can be rewritten as 



k 2 dku) k l{u k - a Q , t) 



f 

Jo 



u! 3 du) 7(a) — a Q , t) . 



(17) 

Eq. (17) has a rather precise physical meaning: it ac- 
counts for the temporal buildup and decay of correla- 
tions between dot states and photons, and thus plays the 
role of a memory function for the reservoir. Assume that 
tw S> 1 which can be easily fulfilled on the time scale 
of the system's dynamics (since a D <~ E a » 8). Then 
7(a) — a> , t) is a function which is strongly peaked around 
a — lu . We thus can approximately replace in Eq. (17) 
a 3 with a 3 . Then, 47 



S (—' 

\ c 



duj 7(a — u> , t) 

n 



dUJJ(UJ — U) ,t) = 7T I 



V c 



) , (18) 



where in the last line we have extended the integral 
to —00. Eq. (18) is the expression we were seeking 
for. It shows that generalized scattering rates can be 
computed without invoking the adiabatic approximation 
which would be problematic in view of the choice of 
asymptotic states. Rather it suffices for the problem of 
our concern to assume that the reservoir memory is suffi- 
ciently short-lived to be approximately treated as a delta 
correlation on the time scale of the system dynamics. 
For the ^-integrals in Eq. (16) we assume 



f 47r 
/ dn(e* ka e s )(e ka e' s ) = (1-/3)— S s 
Jngn„ 6 



nen, 



dQ(e* ka e s )(e k(7 e' s ) = /?y (e*e s )(e ff e' s ) ,(19) 



where two approximations have been employed: first, 
for the unmonitored decay we have extended the inte- 
gral over the complete sphere and have accounted for the 



missing segment fl Q through a reduction factor 1 — 0. 
Accordingly, we have introduced a factor for the moni- 
tored segment Q a , and have assumed that within o the 
polarization vectors are approximately e ka = e CT , where 
e± cx ei ± ie.2- 

With these expressions we can cast the equation of 
motion for the reduced density-operator p d (henceforth 
we will omit the superscript d) to a master equation of 
Lindblad form, 45 ' 53 



-i(H oS p- 



(20) 



where the Lindblad operators Li for the various scat- 
tering channels are listed in table I and H c s = H d — 
i; L\Li is a non-Hcrmitian Hamiltonian accounting 
for the dot states and out-scatterings, which can be sim- 
plified to H cS = B d Q - i£(2 |3)(3| + £ s \s)(s\). 



E. Unraveling of the master equation 

The decomposition of the master equation (20) into 
out-scatterings, H e s, and in-scatterings, second term 
on the right-hand side, allows for a convenient solution 
scheme through unraveling: 35,53 suppose that the system 
is initially in state p(°). Then, p[ a) = U eS (t)p^U^ s (t) 
accounts for the time evolution subject to the condi- 
tion that no scattering has occurred, with U c s(t) = 
cxp —iH e fft. Apparently, this term alone does not pre- 



serve the norm (t) = trp^ since it only describes the 



TABLE I: Lindblad operators for the unmonitored and mon- 
itored decay channels. We use k 2 = (1 — 0)Y and k' 2 = /3L, 
with L the Wigner-Weisskopf decay rate given in Appendix 
A. 



decay channel 


unmonitored decay 


monitored decay 






(polarization 0) 


biexciton 


L :i ^s = k s)(3 


L^ = k»<3| 


exciton 


L s ^o = k|0>(s| 


L% = K'\Q){a\ 
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situation where the system remains unscattered. For un- 
monitored scatterings this lack of norm conservation can 
be restored by introducing a further term which accounts 
for in-scattcrings. The probability for a scattering 
at time t is proportional to the decrease of p(°'(t), i.e., 
p (1) {t) = d(l - p(°\t))/dt, and we obtain 

e*p(°)(t)„<°> + f dt'p^{t')U{t,t')p^U{t',t). 
Jo 

(21) 

Here, /Z 1 ' = i^p^^ij /tr(.) is the density operator 
after scattering and the denominator ensures trp = 1. 
Because of the incoherent nature of scatterings this 
term becomes a mixture of all possible final states. Fi- 
nally, U(t, t') describes the ensuing system's propagation, 
which may again be solved through unraveling (see, e.g., 
Sec. III). ' 

It is important to realize that the mixture in Eq. (21) 
is due to our ignorance about the state of the system and 
reservoir. As the scattering can occur at any time with 
probability p^(t), we have to sum in Eq. (21) over all 
possible histories. A different situation arises if we ac- 
quire additional information about the state of the sys- 
tem by monitoring the reservoir, e.g., through photo de- 
tection. Here, the detection of a photon at time t allows 
to reduce the density operator from a mixture to a single 
term p — » L i( oij/tr(.), where i corresponds to the out- 
come of the measurement, i.e., according to the photon's 
energy and polarization. Note that the monitored scat- 
terings are described identically to the unmonitored ones 
with the only exception of the aforementioned reduction 
of p due to our increase of information. This similarity 
is not incidentally, and can be understood more deeply 
within the framework of pointer states. 4 



F. Photon correlation statistics 

Next, we show how to compute the two-photon corre- 
lation function g^ from the biexcitonic cascade decay. 
We assume that initially the system is in state 3. g^ is 
then obtained through a series of repeated measurements, 
where only those contribute in which photons from the 
same cascade decay are detected by fl s and fl a . Accord- 
ingly, in the following we can disregard photon emissions 
into unmonitored segments. 



1. Biexciton decay 

With the initial density operator p^ = |3)(3| the time 
evolution conditional to no photon emission can be eval- 
uated to = e~ 2r *|3)(3|. When a photon with po- 
larization a- is detected at time t, the density operator 
accordingly has to change to |+)(+|, i.e. a wavefunction 
purification 20 due to the measurement occurs. 



2. Exciton decay 

Since |+) is not an eigenstatc of H% the system will 
start to oscillate between |±). Again we can unravel the 
corresponding master equation. It turns out to be conve- 
nient to use the Pauli-matrix representation of the circu- 
lar basis, Eq. (3), within which the effective Hamiltonian 
reads H c g = (E a — i^)t + |cr 2 - Hence, 

e -iH e{t t = e -i {Eo -iL )t{cos ft^ft (22) 

accounts for the conditional time evolution of no photon 
emission. The corresponding density operator is of the 
form pf ] = |*(t))<*(t)|, where 

= e - iH - B *\+) = e-^-'S^cos ^ |+)+ sin y |-)) 

(23) 

is a state whose norm decreases with time, || ^(t)) || 2 = 
e~ r *. From Eq. (23) we can draw two conclusions: first, 
the probability distribution to detect a second photon 
with any polarization exhibits a mono-exponential de- 
cay without any coherence beats; on the other hand, 
the probability distribution for the detection of a pho- 
ton with well-defined polarization er + is |(o+|\Ij(i))| 2 = 
e~ r *cos 2 y, which indeed exhibits the characteristic 
quantum-coherence beats. Thus, quantum-coherence 
phenomena arc indeed observable in appropriately de- 
signed incoherent optical experiments. 

G. Photon entanglement 

We are now in the position to analyze the shortcoming 
of Fermi's golden rule in describing the problem of our 
present concern. Let us first discuss the biexciton decay. 
From Fig. lb we observe that for this decay the dot-level 
scheme effectively reduces to a A-type scheme, with 3 the 
upper state and s = 1, 2 the lower, approximately degen- 
erate ones. Because of the final-state ambiguity within 
the photon emission process no 'decision is taken' regard- 
ing the final state but rather a quantum-mechanical su- 
perposition is formed: the underlying state corresponds 
to an entanglement between the emitted photon and the 
exciton. Importantly, for a reservoir with sufficiently 
short-lived memory, a good approximation for the system 
under concern, photon and exciton become decoupled on 
a short timescale, thus allowing for the master-equation 
approach employed. 

When the exciton decays through emission of a second 
photon, the photon-exciton entanglement is transferred 
to the two photons. Provided that the photons prop- 
agate along directions where the optical selection rules 
apply, i.e., sectors Q s and Cl a , the two photons are in an 
entangled state; indeed, in Refs. 37,39,40 the authors pro- 
posed such biexciton decay as an entangled two-photon 
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FIG. 2: (a) Schematic band diagram of the proposed struc- 
ture, (b) Level scheme of the spin-degenerate electron states 
a and the charged-exciton state 3 in the dot. 



tions in the contact. As will be shown in the following, 
the framework developed in the previous section can be 
easily extended to describe such propagation in a solid- 
state environment, and we will provide a comprehensible 
theory accounting for the complete cascade process of: 
the buildup of three-particle coherence through tunnel- 
ing; the swapping of quantum coherence to spin entan- 
glement through dephasing and relaxation in the reser- 
voirs; and finally the process of disentanglement through 
spin-selective scatterings. Since the main emphasis of our 
work is on the identification of the basic schemes under- 
lying the buildup and decay of entanglement, we will rely 
on a simplified description scheme of environment inter- 
actions which will allow us to derive analytic expressions 
throughout. 

A. Level scheme 



source. Finally, since the quantum properties of the ini- 
tial photon-exciton system are directly transferred to the 
two photons it is not necessary to detect the first pho- 
ton before the emission of the second one, 56 as one might 
naively assume from our discussion of the quantum-jump 
approach. Accordingly, the proposed scheme would also 
work for photo detectors placed far away from the dot. 

III. SPIN ENTANGLEMENT 

In our second example we discuss the optically trig- 
gered spin entanglement of electrons. Consider the sys- 
tem depicted in Fig. 2a which consists of a single quan- 
tum dot embedded in a field-effect structure. We assume 
that the dot is initially populated by a single surplus 
electron, which can be achieved by applying an external 
bias voltage such that an electron is transferred from an 
nearby n-type reservoir to the dot; 13,15 further charging 
is prohibited because of the Coulomb blockade. Optical 
excitation of this structure then results in the excita- 
tion of a charged exciton, i.e., a complex consisting of 
two electrons and a single hole; 13-15 appropriate tuning 
of light polarization and frequency allows to selectively 
excite the charged-exciton groundstate, where the two 
electrons have opposite spin orientations. Since within 
the field-effect structure the charged exciton is not a sta- 
ble configuration, in a consequent step one electron and 
hole will tunnel out from the dot to the nearby contacts; 
here, the system can follow two pathways, where either 
the electron in the dot has spin-up and the one in the 
reservoir spin-down orientation or vice versa. Similarly 
to the scheme discussed in the previous section, this final- 
state ambiguity results in a maximal spin entanglement 
between the electrons in the dot and reservoir. 

However, as compared to the biexciton decay addi- 
tional difficulties arise because the tunnel-generated elec- 
tron and hole do not propagate freely, as photons would 
in the corresponding scheme, but are subject to interac- 



Our model system comprises the dot, H^, and the 
reservoir degrees of freedom for electrons and holes, H^. 
The dot is characterized by the spin-degenerate single- 
electron states a with energies e CT , and the charged- 
exciton state 3 with energy 63, i.e., 

H d = £ e»H+ e3 |3)(3|. (24) 

<T=1,2 

The electron states in the reservoir are labeled by the 
continuum of quantum numbers k and the spin index cr, 
where the fermionic field operator c ka creates an electron 
with energy e| ff . Correspondingly, d\ creates a hole in 
state k; note that the quantum states k for electrons and 
holes arc different. Then, 

ka k 

In the following we shall not specify the details of the 
reservoir states k, but rather pursue a simplified de- 
scription scheme which only grasps the essential fea- 
tures. To this end, we introduce the fermionic field 
operators C a (uj) = J2 k Ckcr S(ui — e ka ) for electrons and 
D(ui) = J2 k dk S(uj — e k ) for holes, which average over all 
quantum states k with the same energy uj. The following 
anticommutation relations apply, 

{G ff («),C^(w')} = {Cj(w),c£,(w')} = 
{C„M, Cl(cj')} = - u/) , (26) 

with g^(ui) = J2k e% a ) the electron density of states. 
Analogous relations apply for holes. 

B. Tunneling 

Since within the field-effect structure the charged ex- 
citon is not a stable configuration, in a consequent step 
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one electron and hole will tunnel out from the dot to 
the nearby contacts. For simplicity, we assume combined 
electron-hole tunneling, although our results would not 
be significantly modified in case of separate electron and 
hole tunneling. 44 Then, 



V, 



kk'a 



i fefe ^cL4,|a)(3|+h.c.) , (27) 



where t kk /_ a is the matrix element associated to the tran- 
sition from the charged exciton state 3 to a through 
emission of one electron-hole pair. In Eq. (27) we have 
assumed that the total electron spin is a good quan- 
tum number which is conserved in the tunneling process. 
We furthermore assume that t kk i _ a only depends on the 
electron and hole energies e ka and ejj, , respectively, i.e., 
tkk'M = i{e ka ,€ k ,). Hence, the Hamiltonian of Eq. (27) 
can be rewritten as 



Vt = ~J2 f d^edtoh {t{u: ei uo h )Cl{Lo e )D\Lo h )\a){i\+h.c) 

(28) 

Note that all this simplifying assumptions could be easily 
lifted, however, at the price of more cumbersome nota- 
tion. 



1. Master equation 

Assume that the initial density operator is |3)(3|. 57 For 
low temperatures and sufficiently high tunneling rates 
we can safely neglect phonon processes and radiative de- 
cay in the dot, and tunneling becomes the only relevant 
scattering channel. In our theoretical approach tunnel- 
ing is described, analogously to the photon emission of 
Sec. IIC, within lowest-order time dependent perturba- 
tion theory; in addition, we perform the adiabatic ap- 
proximation 2 and replace j(fl,t) by 7r<5(f2). Since the 
hole enters with a high excess energy into the contact it 
immediately suffers an inelastic scattering, which guar- 
antees that tunneling is an irreversible process. This 
is taken into account by tracing over the hole degrees 
of freedom and neglecting terms tr^p-D^ (w')-D(w) = 0. 
Then, tr hP D(uj)D^(u') = g h (cj)S(uj - w'). Within this 
framework we arrive after some straightforward calcula- 
tion at the master equation 



pt ~ -i(H eSPt - p t Hl s ) dujduJ ' 

))ct(u;)\a)(3\p t \3)(a'\CA^), 

(29) 

which includes tunneling as the only scattering channel. 
Here, T a (u>,uj') = J g h (uj)duj t(u>, uj)t{uj' , uj)w5(uj + lu + 
e s — £3) and we have assumed that t is real. Furthermore, 



H cti = H - i^ a Jdu>du>'f a (u,u')\3)(3\C^u')Ci(u) is 
an effective Hamiltonian which accounts for the dot and 
reservoir states, H a — + H^, and outscatterings. 

Next, we compute the term H c ffp. We assume that 
the electrons in the reservoir can be described as an elec- 
tron gas at zero temperature with total wavefunction 
rifecr ®( e F — e ka ) c k<r\0), where Heaviside's step function 9 
ensures that only states with an energy below the Fermi 
energy ep are populated. Hence, 



CM)CI{lo)p 

= 9{u> - e F ) (&{«,)6{w - u') - Ct^C^J)) p 
= 6{u-e F )5(Lu-uj')g e a {uj)p. (30) 

The term Cl(oj)C a (ui')p in Eq. (30) accounts for a pro- 
cess where an electron is exchanged between reservoir 
and dot. In the following we will neglect such processes 
since they only contribute to an energy renormalization 
but not to scattering and dephasing. From Fig. 2 we 
infer that the energy of the tunneled electron must be 
larger than ep and smaller than some cutoff energy u> c , 
which is determined by the kinematics of the tunneling 
process. To simplify the following calculation, we assume 
that within the energy window ep < lu,lu' < lu c the quan- 
tities T a (oj,uj') = T and g%{us) = g% are approximately 
constant. Then, H eS = H Q - iT' t \?>){2>\ with 



dl{u)T = 5o A T = y , (31) 



where T' t denotes the single-electron tunneling rate and 
A = u c - e F . 



2. Unraveling of the master equation 

Within this scheme we can solve the master equation 
subject to the initial density operator |3)(3| through un- 
raveling, 



pt S e- rt |3)(3| +T t ( dt'e- r ^'U{t,t')p^U{t',t) , 
Jo 

(32) 

where T t = 2T' t is the charged-exciton tunneling rate and 
U(t, t') the time evolution operator accounting for the 
propagation of the tunneled electron (see below). The 
density operator after tunneling is p^ = \^>i){^i\, with 



|* 1 ) = (2^A)"^ rcl{u)\a). 



(33) 



Eq. (33) is an important and non-trivial result. First, it 
demonstrates that despite the incoherent nature of tun- 
neling and hole relaxation the electron system can be 



10 



TABLE II: Contributions to the density operator corresponding to the different steps of the cascade process, i.e., for the initial 
charged exciton, the system after tunneling, after dephasing, and after disentanglement. The last column shows the time 
evolution of the corresponding probabilities. 



step 


density operator 


time evolution of probability 


before tunneling 


|3)<3| 




p<°> = -r tP (°) 


after tunneling 


(2 ff sA)- i x;„/ r; dwdw'cjM 


\*)(a'\C a (u') 


pd) = r tP P» - i>w 


after dephasing 


(2g e A)- 1 Z^ Qdwdoj'CKuj) 


\9){9'\C a (u')S(u-u') 


p (2) = r dP « - 2 K 2 r>< 2 > 


after disentanglement 


(2 ff SA)- 1 E„.C^'Cj(w) 


\a){a'\C*(u')6(u-u')S„, 


p< 3 > = 2k 2 iV 2) 



described in terms of wavefunctions; note that the de- 
tection of the hole would even allow to purify this wave- 
function, 20 which might be of relevance when initially p 
is not equal to |3)(3|. Secondly, a closer inspection of 
Eq. (33) reveals that the spin part C + \— ) + C_|+) is a 
maximally entangled state of the electrons in the dot and 
reservoir. We emphasize that this maximal entanglement 
is independent of the spin basis, which guarantees that 
our scheme is not deteriorated by possible polarization 
anisotropies of the dot states, e.g., fine-structure split- 
tings. 

C. Dephasing 

After tunneling the electron in the reservoir propa- 
gates in presence of scatterings, as described by U(t, t') 
in Eq. (32). Quite generally, we assume that the or- 
bital degrees of the reservoir electron are subject to much 
stronger interaction channels, e.g., phonons, than the 
spin degrees, as evidenced by the long measured spin 
lifetimes (~ns) in n-doped semiconductors. 58 For that 
reason, let us first consider an elastic electron scattering 
which does not depend on spin, i.e., Lindblad operators 
of the form 

LM = r | E ^™ (34) 

with T d the scattering rate. Hence, the effective Hamil- 
tonian reads H c s = H a — | J du; (lo)L(lo). With 
/ dwL(w)|\l/i)(\E , i|L t (cj)/tr(.) the density operator after 
scattering and L(w)|*i) = (2^A)"5 £ ff C„(u)\d), the 
master equation accounting for elastic spin-unselective 
scatterings can be solved through unraveling. Suppose 
that the tunneling has occurred at time t', correspond- 
ing to one of the possible histories in Eq. (32). Then, 

UtyJ^Utfj) = e- r ^*- t ')e- iif °*|* 1 )(*i|e iff <' t 

+Y d / die-r^-^Uit^pfuilt) (35) 
Jt' 

is the corresponding conditional density operator at later 
time. Here, 



pp ) = (2 5 ^A)- 1 ^ duCt{u>)\a){a'\CAu) (36) 

(7(7 

is the density operator after the elastic scattering and 
U(t, t) accounts for the propagation of the scattered elec- 
tron. In comparison to Eq. (33) the density operator of 
Eq. (36) is diagonal in cj, i.e., the elastic scattering has 
led to a destruction of the phase coherence (dephasing). 
However, the spin part still shows the same degree of 
entanglement, where similar conclusions would apply for 
inelastic but spin-independent scatterings. Thus, the de- 
cay of an optically excited charged-exciton indeed gener- 
ates a robust spin entanglement between the electrons in 
the dot and reservoir. 



D. Disentanglement 

We finally comment on the process of disentanglement. 
In fact, any scattering channel which couples with un- 
equal strength to the spins or affects only one spin ori- 
entation is responsible for such entanglement decay. To 
simplify our following analysis, we introduce operators 
acting in the reservoir subspace 

(^A)- 1 r dwCjMCVM - V aa , , (37) 

J ep 

which are averaged over the orbital degrees of freedom 
uj. Without specifying the details of such spin-selective 
scattering, we introduce the generic Lindblad operator 

L = 17 ( (1 + k)Vh + (1 - k)T 2 2 ) , (38) 

where T s is a scattering rate and n a factor determining 
the asymmetry of the coupling to spin-up and spin-down 
electrons. As briefly sketched in Appendix B, the corre- 
sponding master equation can be solved analytically and 
we arrive at 

U(tX) P {" ) U(t\t)^e-^ t - t '^'P^)^\ 

<J<j' 

+ (1 - e- 2K2T ^-^)J2v a<7 \a){a\ . (39) 
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FIG. 3: Schematic representation of the various contributions 
p' l '(t) of Eq. (21) as discussed in the text. Times are mea- 
sured in units of the inverse tunneling rate Fj -1 ; we assume 
r d = 2Y t for the spin-unselective and 2n 2 T a = 3x 1CT 3 Y t 
for the spin-selective scatterings. The black dashed line shows 
the sum of p^(t) and p^it) which corresponds to the degree 
of spin entanglement. 

Here, the first term on the right-hand side accounts 
for the unscattered state, whose probability decays with 
time, and the second term for the disentangled state af- 
ter scattering. Because of our assumption for the Lind- 
blad operators, Eq. (38), the density operator after dis- 
entanglement still exhibits a classical correlation where 
the spins of the electrons in the reservoir and dot are 
antiparallel; apparently, it would require additional spin- 
flip processes to completely decouple the two electrons. 
See also Ref. 44 for a discussion of how such spin entan- 
glement could be measured experimentally and possible 
quantum information applications. 

E. Cascade decay 

We finally list in table II the various contributions of 
the cascade decay to the density operator, i.e., corre- 
sponding to the step before tunneling, after tunneling, 
after dephasing, and after disentanglement. The prob- 
abilities for the corresponding contributions are deter- 
mined by our unraveling procedure and can be obtained 
as follows: initially, the system is in the charged-exciton 
state, i.e., p(°\0) = 1; because of tunneling decays 
mono-exponcntially; correspondingly, the increase of the 
contribution after tunneling p^ is proportional to the 
decrease of p^ , which again diminishes because of de- 
phasing, and so forth. Within this scheme, we can easily 
obtain the time evolution of the various probabilities as 
shown in the last column of table II, with the initial con- 
ditions p(°)(0) = 1 and pW(0) = for i = 1,2,3. Fig. 3 
sketches the various contributions for representative val- 
ues of Td and T s . 

IV. CONCLUSION 

In conclusion, we have analyzed single scatterings in 
single semiconductor quantum dots for two representa- 



tive examples. First, for the cascade decay of a biexciton 
it has been shown that an exciton fine-structure splitting 
results for appropriately chosen polarization filters in an 
oscillatory behavior of the two-photon correlations. In 
our second example we have investigated the tunneling 
decay of a charged exciton inside a field-effect structure 
into a photo current, and have shown that the spins of the 
electrons in the dot and reservoir become entangled. We 
have discussed that this entanglement is robust against 
dephasing and spin-unselective scatterings, and thus ben- 
efits from the long spin lifetimes in semiconductors. 

Our theoretical approach has been based on a master- 
equation approach in Lindblad form in conjunction with 
the quantum-jump approach for the description of pho- 
ton measurements. We have seen that the description 
scheme in terms of Lindblad operators, suitable for the 
environment couplings of our present concern, is supe- 
rior over related quantum-transport descriptions based 
on the Markov and adiabatic approximations. 2 First, no 
assumption about the final states of a 'scattering' has to 
be invoked, and it suffices to assume a sufficiently short- 
lived memory kernel of the reservoirs instead. Secondly, 
Lindblad operators associated to measurement, e.g., pho- 
ton detection, allow for the application of the quantum- 
jump approach and provide a means to describe single- 
system dynamics. Finally, the distinction of out- and 
in-scatterings in the master-equation makes possible a 
flexible solution scheme through unraveling, which ap- 
pears to be highly elegant in particular for the cascade 
processes discussed in this work. We expect our find- 
ings to be useful for the analysis of time-resolved single- 
system measurements and for the simulation of quantum- 
information applications. 
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APPENDIX A: WIGNER-WEISSKOPF DECAY 
RATE 

In this appendix we show how to evaluate contribu- 
tions of the form 27T Y^ka- 9k,os9k,crs< S(uj k - w ). First, 
the polarization vector of a photon with wavevector k 
and polarization a can be expressed as 47 

e>fc !+ = — -= (cos 9 cos (f> — i sin </>, cos 9 sin <j> + i cos <f>, — sin 
v2 

&k- = — j= (cos 9 cos 4> + i sin <f>, cos 9 sin <j) — i cos <j), sin 



The integral over the angular-dependent part can be eval- 
uated to 
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sm9d9\e* ka e s \ 2 



flit f-~K 

/ # / I 

/ d<j> sin^(eLei)(e fcCT e 2 ) -0, (A2) 
Jo Jo 



i.e., only terms with s = s' give a nonzero contribution. 
Finally, 

T = 2ir^2 \9ka\ 2 5{uj k - ui ) 

= 2,f,2, r 3^(^)' = ±gi,(A 3 , 

is the Wigner-Weisskopf decay rate for a dipole radiator 
embedded in a medium with refractive index n. 



with the vector cr = (ci, 02, (T3). From trp = 1 we find 
Mo = 5 an d P — P implies that m is real. Analogously, 
a generic Lindblad operator is of the form L = a^t + 
a • cr, where a = a + ia^ and a = a' + ia" can be 
decomposed into real and imaginary parts. The master 
equation corresponding to this Lindblad operator is of 
the form 



= \([Lp,tf] + [L,ptf}) , 



(B2) 



which can be easily evaluated by using the usual commu- 
tation relations for Pauli matrices. We arrive after some 
lengthy calculation at 



APPENDIX B: MASTER EQUATION FOR 
TWO-LEVEL SYSTEM 

In this appendix we discuss the master equation of a 
generic two-level system. Since the Pauli matrices <7j 
together with the unit matrix 1 form a complete basis, 
we can expand the density operator as 



p = uq! + u ■ cr . 



(Bl) 



ii = 2( (a' x a") - (a' a" - a V) x u 

-\a\ 2 u + (a'u) a' + (a"u) a"^j . (B3) 



Then, for the Lindblad operator L = (1 + na 3 ) of 
Eq. (38) for disentanglement we find W1.2 = — 2n 2 Tui t 2 
and u 3 = 0, which can be easily solved to arrive at the 
final result of Eq. (39). 
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